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measured aircraft position vector 
time from start of simulation run 
smooth-commanded aircraft velocity vector 
total commanded aircraft acceleration vector 
measured aircraft velocity vector 
commanded angle of attack 
commanded vertical flightpath angle 
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<j> c commanded roll angle 

commanded horizontal flightpath angle 
co measured aircraft angular velocity vector 

co c smooth-commanded aircraft angular velocity vector 

*c total commanded aircraft angular acceleration vector 
w open-loop commanded aircraft angular acceleration vector 

d) cc closed-loop commanded aircraft angular acceleration vector 
Acronym 

TAFCOS total aircraft flight-control system 
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APPLICATION OF THE CONCEPT OF DYNAMIC TRIM CONTROL 


TO AUTOMATIC LANDING OF CARRIER AIRCRAFT 
G. Allan Smith and George Meyer 
Ames Research Center 


SUMMARY 


The availability of the airborne digital computer has made possible the 
practical implementation of the concept of a balanced feed-forward and feed- 
back automatic aircraft f lightpath-control system. The concept is explained 
and simulation results are presented. The essential nonlinear force and moment 
characteristics of the aircraft are collected in tabular form as dynamic trim 
maps and inverted to give a feed-forward command signal path which, in series 
with the actual aircraft, provides essentially an identity transfer function. 
With perfect modeling and no disturbances this would provide perfect trajec- 
tory control. Feedback loops are closed around this linear path to compensate 
for disturbances and imperfect modeling. Simulation results and a flight test 
have shown that only a small fraction of the total drive signal is required 
from the feedback while the major portion is provided by the feed-forward 
control. 

This concept has several advantages. It allows a straightforward design 
for aircraft with highly nonlinear characteristics since conventional linear 
techniques can be employed to design the feedback loops around the linearized 
feed-forward path. The use of the complete aircraft characteristics in the 
feed-forward path allows an integrated design so that total aerodynamic and 
propulsive forces and moments are combined in a natural way and aid, rather 
than partially offset, each other as often occurs in conventional designs. 

The explicit aircraft aerodynamic characteristics available in the feed-forward 
signal path allow direct limiting of commanded angle of attack and various 
accelerations and rates so that a smooth executable trajectory is commanded 
regardless of excessive or inadvertent commands from an air-traffic-control 
system or other input. Required digital computer speed and memory is substan- 
tially less than that required for conventional designs of comparable capa- 
bility. 

Simulation results are presented for an application to automatic landing 
on an aircraft carrier. 


INTRODUCTION 


Precise flightpath control is required over the approach trajectory to 
land an aircraft on an aircraft carrier in heavy seas. The touchdown point is 
restricted to a small area on a deck that is pitching and heaving due to wave 
action. Furthermore, the flightpath is continually disturbed by severe atmo- 
spheric turbulence generated in part by the carrier superstructure. Landing 



the aircraft is particularly difficult at night and is generally acknowledged 
to be the most difficult piloting task encountered in routine operations. 

In an attempt to improve landing performance, the Navy has developed an 
automatic carrier-landing system. The system is based on a conventional auto- 
pilot that is controlled by attitude signals sent by radio link from the car- 
rier and generated from path errors determined by a carrier-based radar. 
Although the system is operational for a few types of aircraft, it has perfor- 
mance deficiencies in strong turbulence and significant touchdown point dis- 
persion. The alternative automatic-control-system design concept presented in 
this report has shown potential for improving landing performance. 

This new design concept was developed at Ames Research Center over the 
past 4 years. The theoretical basis for the system is presented in refer- 
ence 1. This report will explain the system structure developed to apply the 
concept to the carrier-landing problem and present results of simulation 
studies of an automatic carrier-landing system based on the new principles. 

The concept provides for control of the aircraft trajectory directly by 
commanding the aerodynamic and propulsive forces required rather than by a con- 
ventional autopilot that commands aircraft attitude to minimize trajectory 
errors. In addition, the dynamic trim concept allows a single configuration 
to function over the complete flight envelope without the gain scheduling or 
configuration switching that is used in conventional systems based on pertur- 
bations about a number of different static trim points. The concept is termed 
a M total aircraft flight-control system," or TAFCOS. TAFCOS is essentially an 
open-loop, feed-forward system that commands the proper instantaneous thrust, 
angle of attack, and roll angle dynamic trim conditions to achieve the forces 
required to follow the desired trajectory. The dynamic trim conditions are 
determined by an inversion of the aircraft nonlinear force characteristics to 
calculate the required angle of attack and throttle setting, while the com- 
manded roll angle is determined from the required normal and lateral force 
coefficients. The angle of attack is combined with roll angle and flightpath 
angles to give the commanded attitude. This feed-forward design is completed 
by inversion of the aircraft nonlinear moment characteristics to determine the 
control surface angles. 

The concept of open-loop control for a completely calibrated system can 
provide any desired response that is within the physical capability of the 
system. On the other hand, the value of feedback to deal with uncalibrated 
systems or systems with uncertainties and disturbances is well established. 
Therefore the basic open-loop design of TAFCOS is supplemented by feedback 
loops that compensate for lack of knowledge of the aircraft characteristics 
and for imperfect wind estimation, and provide control in the presence of dis- 
turbances. TAFCOS is a balanced system that employs an open-loop controller 
with as much a priori information as feasible about the aircraft force and 
moment characteristics; it uses feedback to provide tight control in the face 
of limited knowledge of plant characteristics, measurement inaccuracies, and 
disturbance inputs. 

Advances in avionics technology have produced airborne digital computers 
of the type needed to carry out the calculations required for TAFCOS. 
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Extensive two-dimensional tables of the aircraft force and moment characteris- 
tics are used for the inversion calculations, and numerous matrix multiplica- 
tions are required to transform the force and moment vectors between various 
coordinate systems. Computational feasibility of the concept has recently 
been demonstrated by a successful flight test of TAFCOS in a STOL aircraft 
(ref. 2). A Sperry 1819A airborne digital computer was used. Only a small 
portion of the computer capacity was required. 

During the past 3 years several simulations of the current operational 
Navy automatic carrier-landing system have been carried out on the piloted 
simulators at Ames Research Center. They were conducted jointly by Ames 
Research Center, the Naval Air Test Center, and the Naval Air Development 
Center to investigate pilot acceptance of various system modifications de- 
signed to improve touchdown point dispersion under turbulent conditions 
(ref. 3). The same computer simulation model of the A-7E aircraft used by the 
Navy for these simulations was used to represent the aircraft in the TAFCOS 
simulations that are the subject of this report. The model includes the non- 
linear aerodynamic force and moment equations and nonlinear effects of the 
engine, including hysteresis in the throttle actuator. 

The simulations presented here were carried out to demonstrate that the 
TAFCOS concept provides substantially improved performance for the carrier 
landing task. The results are presented in four sections to show the critical 
responses that must be examined to evaluate automatic control-system perfor- 
mance. The first results show all the important variables for a tip-over and 
landing maneuver without atmospheric turbulence. This allows the attitude and 
thrust control interaction to be clearly observed and establishes a perfor- 
mance reference for the task of tracking the moving deck. The next results 
show the performance for a landing with atmospheric turbulence. 

The quality of an automatic control system is established by both its 
ability to regulate to a reference flightpath in the presence of input distur- 
bances and to respond to maneuver commands. Therefore, the final two sections 
of results show the system response to a series of severe gust inputs and 
response to a sequence of maneuver commands that include altitude changes and 
lateral flightpath control. 


AUTOMATIC CARRIER-LANDING TASK 


The automatic carrier-landing situation is illustrated in figure 1. The 
carrier normally heads into the wind at a speed such that the relative wind is 
parallel to the deck at about 30 knots. The aircraft is first directed to a 
marshall point or to a holding pattern. Individual aircraft are then cleared 
to the approach course at 2-min intervals. For the A-7E aircraft considered 
in this simulation study, the aircraft is commanded to hold a steady velocity, 
with respect to the wind of 129 knots, which is 99 knots with respect to the 
carrier. 


For this simulation a level-approach course is commanded at an elevation 
of 152.4 m (500 ft) above the carrier-deck reference until a 3.5° glide slope 
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to the ideal touchdown point is intercepted. The ideal touchdown point is 
located 70 m (230 ft) aft of the carrier center of pitch. The aircraft must 
touch down within ±18.3 m (±60 ft) longitudinally and ±4 . 6 m (±15 ft) laterally 
of that point. 

The carrier is equipped with a tracking radar that determines the air- 
craft position with respect to the carrier. A shipboard inertial platform 
measures the instantaneous elevation of the ideal touchdown point due to deck 
pitch and heave with respect to the undisturbed carrier deck reference posi- 
tion. As shown in figure 1, after tip-over, the aircraft is commanded to 
descend on the glide slope until 12 sec before touchdown. At that time the 
aircraft is commanded to follow the instantaneous displaced position of the 
touchdown point. 

An approach that misses the required touchdown area or an arresting wire 
is termed a bolter and the aircraft is commanded to follow a recovery path back 
to the reference altitude of 152.4 m (500 ft) before making another attempt. 

The most critical portion of the landing sequence is the final 12 sec, 
during which time the aircraft is commanded to follow the instantaneous verti- 
cal position of the touchdown point. For this simulation, the deck motion was 
taken as the sum of two sinusoidal displacements due to heave and pitch, the 
former at a frequency of 0.1 Hz and an amplitude of ±1.2 m (±4 ft), and the 
latter at a frequency of 0.096 Hz and an amplitude of ±1°, which is ±1.2 m 
(±4 ft) vertical displacement at the 70.1-m (230-ft) distance from the center 
of pitch. This can produce a maximum touchdown point excursion of ±2.4 m 
(±8 ft). Also, during the final portion of the run, strong atmospheric turbu- 
lence, called burble, is encountered in the carrier wake. This turbulence is 
caused by the interaction of the carrier superstructure with the wind; it was 
modeled as a combination of vertical and horizontal gust patterns in time and 
space as a function of carrier deck attitude plus random superimposed turbu- 
lence, Peak total turbulence was about ±1.8 m/sec (±6 ft/sec). 

The primary concern in this simulation was to investigate aircraft 
landing-point dispersion due to the effects of pitch and heave deck motion 
and atmospheric turbulence. Therefore, 160 landing simulation runs were 
made with various random turbulence patterns and with different deck motion 
component phases at touchdown. 

It should be noted that a carrier landing is accomplished by flying a 
direct path to the deck without any landing flare. 

For the currently operational Navy automatic carrier-landing system the 
aircraft is under manual control during a holding pattern, or after a missed 
approach during wave-off. However, TAFCOS has the potential for automatic 
control during such situations, so simulation runs were made of additional tra- 
jectories that incorporated lateral maneuvers, altitude changes, and severe 
gusts. These runs included maneuver commands and disturbances that illustrate 
the system response over a much wider flight envelope than the conventional 
carrier approach path. 
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An important part of any automatic trajectory-control system is the 
instrumentation required to measure the aircraft position with respect to the 
commanded trajectory. The current Navy system uses a ship-mounted tracking 
radar in conjunction with a beacon transponder on the aircraft. Such a system 
was therefore simulated for TAFCOS; it included the usual loss of radar infor- 
mation when the aircraft is within about 150 m (492 ft) of touchdown. Of 
course, TAFCOS receives the radar position information by data link from the 
carrier and processes it in the airborne digital computer. This is in con- 
trast to the operational system, which processes the radar data in a shipborne 
computer and sends aircraft attitude commands to the approaching vehicle. The 
radar position and the velocity and acceleration data used by TAFCOS were cor- 
rupted by appropriate noise inputs. 


AUTOMATIC CONTROL SYSTEM 


The command structure of TAFCOS is entirely different from that of a con- 
ventional control system. It employs dynamic trim maps which represent an 
inversion of the nonlinear six-degree-of-freedom equations of the aircraft. 

The trim maps take the form of extensive two-dimensional tables of force and 
moment data. They represent a calibration of the aircraft so that trajectory 
commands can be processed in an open-loop fashion to yield surface and thrust 
controls. Closed-loop state feedback is employed to balance trim map inaccu- 
racies and measurement uncertainties. 

Figures 2 and 3 will be used to describe the configuration and functions 
of TAFCOS. They illustrate the signal flow paths that will be traced tc show 
what operations are carried out by the system. The details of how these opera- 
tions are accomplished are presented in appendix A. Many of the quantities in 
figures 2 and 3 are vectors which, for calculation, are resolved into compo- 
nents in appropriate coordinate systems. For example, the trajectory is pre- 
sented in an inertial frame, the major aerodynamic forces act in a frame 
determined by the aircraft velocity with respect to the air mass, and the air- 
craft control surfaces and engine thrust are expressed in an aircraft body 
reference frame. Therefore, a characteristic of TAFCOS is the repeated reso- 
lution of forces, moments, accelerations, velocities, and positions from one 
reference frame to another using matrix transformation methods. Figure 2 
illustrates only the open-loop elements and does not include the closed-loop 
feedback employed in the complete system. Figure 2 will be used for a brief 
explanation of the fundamental concept of TAFCOS without consideration of many 
secondary, although important, details that are treated more completely in the 
discussion of figure 3 and in appendix A. 

The desired trajectory in figure 2 enters the system as a series of 
straight, circular, or helical arc geometric segments described by their 
length, orientation, and airspeed. Angles of climb or descent for all seg- 
ments, radii, and arc lengths for circular and helical sections are given. 

From this geometric-data input the trajectory command sequencer produces a 
consistent commanded dynamic sequence of trajectory position, velocity, and 
acceleration vectors in inertial space coordinates. These dynamic commands 
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are rough in the sense that they may require instantaneous changes in velocity, 
direction, and, sometimes, position at the segment junctions. The trajectory 
is defined in space relative to the carrier and at a commanded airspeed. 

The trajectory command generator smooths and limits these rough commands 
to' be consistent with the desired aircraft response. The trajectory command 
generator also performs a coordinate system rotation so that the desired tra- 
jectory commands in inertial space axes are transformed to relative wind veloc- 
ity axes. The output vector V c is a smooth executable commanded accelera- 
tion. Smooth-commanded velocity components are combined with estimated steady 
wind to calculate commanded flightpath angles, T c and ip c . 

In the coefficient development section, V c is multiplied by aircraft 
mass to give the force required to carry out the desired trajectory. The force 
is then divided by dynamic pressure and wing area to give the total lift and 
drag coefficients. These total commanded coefficients, which include engine 
thrust forces as well as aerodynamic forces, are sent to the force trim map. 
Furthermore, the ratio of the lateral to the normal commanded force coeffi- 
cient gives the tangent of the commanded roll angle <j> c . Figure 4 is a graphi- 
cal representation of the force trim map which is seen to be a series of lift 
drag polars for the aircraft, shifted by the effect of engine thrust as indi- 
cated by the thrust coefficient C x c • Inputs to the trim map are the commanded 
lift and drag coefficients Cl c and Cd c * The trim map outputs are the corre- 
sponding commanded angle of attack a c and thrust coefficient C-j. calcula- 
ted by the digital computer, which carries out an interpolation of the trim-map 
data stored in tabular form. This is a dynamic trim map in the sense that 
angle of attack and thrust are calculated for any vector acceleration, not 
just for a conventional trimmed flightpath. A very convenient feature of the 
trim map is the ease with which it allows commands to be limited to selected 
values of a c and C-j c . For this simulation, commanded angle of attack was 
limited to between -2° and +16° and commanded thrust coefficient to between 
0.01 and 0.7. This permits commands over almost the entire flight envelope 
but prevents commands that approach stall conditions. The effect of this 
limiting is clearly seen in figure 8(e) of the results. 

In figure 2, one output of the trim map, commanded thrust coefficient, 

C’j , is converted to commanded thrust and sent to an engine thrust table where 
the digital computer performs an inverse interpolation to calculate the com- 
manded throttle angle. The other output of the trim map, commanded angle of 
attack a c , is combined with commanded roll angle and horizontal and vertical 
flightpath angles to give the commanded aircraft attitude. 

A procedure generally similar to that followed in the trajectory portion 
of the system is then employed in the attitude control section. The attitude 
command generator provides a smooth executable angular acceleration vector com- 
mand d) c . This angular acceleration vector is multiplied by the aircraft 
inertia matrix to give the torques required to carry out the desired attitude 
variations. These torques are then expressed in terms of moment coefficients. 
The moment trim map in figure 2 inverts the aircraft torque equations to yield 
the control surface angles needed to generate the commanded angular accelera- 
tion vector. 
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When the basic feed-forward control of figure 2 is implemented with rea- 
sonably accurate trim maps, the aircraft follows the commanded acceleration 
V c quite well. The output response is essentially linearly proportional to 
the input command. Thus, it can be considered to be a linear acceleration 
command system and linear theory can be used to select the feedback gains used 
for the actual complete control system (ref. 4). 

The basic open-loop concept of figure 2 is augmented in the complete 
TAFCOS configuration by the addition of feedback loops as shown in figure 3. 
Figure 3 shows the underlying structure of figure 2 with additional symbolism 
as well as the closed-loop position and attitude feedback. Only the new addi- 
tions and certain details previously ignored will be discussed in connection 
with figure 3. 

The commanded acceleration V c in figure 2 includes a component to bal- 
ance the acceleration of gravity so it is more accurately a specific force or 
the total force per unit mass required to follow the trajectory under the 
influence of gravity. The symbol f tc is therefore introduced to indicate 
commanded specific force in figure 3. 

In the trajectory regulator of figure 3 the smooth commanded position R c 
is compared with the measured position R and the smooth commanded velocity 
V c is compared with the measured velocity V. The resulting error signals 
are combined with appropriate gains and limits to give a closed-loop specific 
force command f cc . This perturbation command from the trajectory regulator 
is added to the open-loop command f QC to form the total specific force com- 
mand ftc* For most maneuvers, f cc is only about 20% of f QC . This means 
that the major burden of control is borne by the open loop. 

The trajectory command generator and the trajectory regulator are the 
principal parts of the system where gain and limit adjustments are made to 
achieve satisfactory performance, which is characterized by smooth and accurate 
response to trajectory commands, properly coordinated and suitably limited 
control-surface and engine-thrust responses to disturbance inputs, and well- 
executed capture of the initial trajectory and transition between trajectory 
segments. 

As previously noted, the force trim map for this application expresses 
the lift and drag coefficients as functions of angle of attack and thrust. 

For a different application, a three-dimensional trim map was developed with 
flap angle as a third parameter and a three-parameter interpolation routine was 
used in the airborne digital computer (ref. 5). For attitude commands, TAFCOS 
manipulates angles by constructing the corresponding direction cosine matrices 
that represent the angular rotation about the appropriate axis (ref. 6). The 
required angle of attack a c and roll angle are combined with commanded 

flightpath angles F c and by multiplication of their direction cosine 

matrices to give the rough commanded aircraft attitude matrix A aSrc which 
serves as the input to the attitude command generator. The attitude command 
generator provides a smooth commanded attitude matrix A aSc , a smooth commanded 
angular velocity vector u) c , and a smooth commanded open-loop angular accelera- 
tion vector u) oc . The attitude regulator compares the commanded and measured 
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attitudes and angular velocities to form error signals which are combined with 
suitable gains and limits to give a closed-loop perturbation angular accelera- 
tion command £> cc . This is added to the open-loop command u) oc to give the 
total angular acceleration command vector a> c . 

The moment trim map in figure 3 is functionally similar to the force trim 
map. It accepts commanded moment coefficients as inputs and provides the com- 
manded control surface positions necessary to achieve the corresponding angular 
accelerations. These commanded surface deflections are the final output of the 
airborne digital computer. They are sent to conventional control-surface 
servos which produce the actual deflections as influenced by servo lags and by 
surface rate and deflection limits. Additional details of the TAFCOS structure 
are discussed in appendix A in connection with figure 10 which is a further 
expansion of figure 3. 

In summary, TAFCOS is an aerodynamic ally oriented computational flight 
control system that fully incorporates the limitations and nonlinearities of 
the force and moment characteristics, including the propulsion system, so that 
it is possible to maintain desired stall margins and limits on attitude 
response rates and control deflections, and yet exploit the aircraft capabili- 
ties over its full practical flight envelope. 


SIMULATION RESULTS 


Simulation runs were carried- out on a general purpose digital computer 
(IBM 360/67). The results appear as time-history plots of significant vari- 
ables as measured during the approach to the carrier. Measurements are taken 
in a right-handed orthogonal coordinate system with origin at the undisturbed 
center of pitch of the carrier. The carrier is on a straight course with a 
30-knot wind over the deck which is aligned with the X axis of the inertial 
space coordinate system as defined in appendix B. The carrier heave and pitch 
motions are the same as those used for the Navy simulations mentioned previ- 
ously; sinusoidal heave and pitch that produce a maximum vertical excursion of 
the touchdown point of about ±2.4 m (±8 ft). 

The results are arranged in four groups. The first group includes fig- 
ures 5(a-h) which show results for a tip-over and landing with carrier deck 
motion but no atmospheric turbulence or other disturbances. This portrays the 
upper bound of performance level, allows a comparison of smooth variables, and 
illustrates theoretical design considerations. 

The second group includes figures 6(a-f) and 7 which show results for a 
tip-over and landing with carrier deck motion following, atmospheric turbulence, 
radar noise, and throttle hysteresis to give a realistic simulation. One hun- 
dred and sixty of these runs with different disturbance statistics were made 
to collect data on touchdown dispersion. Results from a typical run are pre- 
sented in the figures. 
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The third group includes figures 8(a-f) which show results for a straight 
and level run with no atmospheric turbulence but with a series of sharp, 
severe gusts to show how TAFCOS responds to disturbance inputs. 

The fourth group includes figures 9(a-f) which show results for a modi- 
fied racetrack course with large altitude changes. There is no atmospheric 
turbulence, but a steady 30-knot wind is imposed. These curves show how TAFCOS 
responds to a range of lateral and vertical commands at a commanded airspeed 
of 129 knots for a varying wind orientation to the flightpath. 


Carrier Approach with No Disturbances 

The performance of the A-7E aircraft with the TAFCOS automatic carrier 
landing system for a simplified standard touchdown run with no atmospheric tur- 
bulence is shown in figures 5(a-h). For figures 5(a-e), no hysteresis is 
included in the throttle actuator, and radar-range information is used all 
the way to touchdown. Moreover, the phase of the carrier sinusoidal deck 
motion is such that the touchdown occurs at the peak of touchdown point dis- 
placement from its undisturbed position. These runs thus establish the best 
system performance; they will be compared with other runs that include more 
realistic operating conditions. The aircraft is moving in the direction of 
the positive X and Z axes toward the carrier at the origin. Thus the air- 
craft position coordinates are generally negative in both the X and Z axes. 

Figure 5(a) consists of time histories of the trajectory variables in the 
vertical channel. Tip-over occurs at 8 sec where the actual vertical position 
and velocity show a very smooth transition. To provide anticipation, the rough 
position and velocity commands jump from the level-flight trajectory to a glide 
slope of 3.5°. A section of the initial vertical-position curves, plotted to 
expanded time and altitude scales, is presented to show how effectively the 
trajectory command generator smooths the rough-position command. The smooth 
position, velocity, and acceleration commands are very closely followed by the 
actual aircraft response. It will be noted that although no rough-commanded 
vertical acceleration is provided at tip-over, the trajectory command generator 
does provide a smooth tip-over acceleration command. It is interesting to 
note that the actual acceleration closely follows the smooth command, except 
for a slight initial acceleration at tip-over that is in the opposite direc- 
tion; and is characteristic of the nonminimum phase, right-half plane zeroes 
of the aircraft transfer function. 

Figure 5(b) shows the carrier deck motion and the deck motion following 
commands that are superimposed on the rough-commanded glide-slope trajectory 
during the final 12 sec before touchdown. Because these quantities are mea- 
sured by the ship inertial platform they can be transmitted by radio link to 
the aircraft until touchdown. It will be noted that these signals are atten- 
uated during the first 2 sec so that they build up slowly to their final val- 
ues and do not jump abruptly. It will also be observed that touchdown occurs 
at the peak of carrier motion. One of the most important curves in fig- 
ure 5(c) is the throttle (power lever angle) which shows, except at tip-over, 
a very smooth response consistent with the acceleration requirements of the 
aircraft longitudinal axis. The other variables of figure 5(c) show 
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satisfactory transient response at tip-over and recovery to steady values 
before taking on sinusoidal variations during the deck motion following just 
before touchdown. For this run the aircraft passed over the carrier ramp 
with a clearance of 5.2 m (17.1 ft) and landed 0.61 cm (0.02 ft) beyond the 
ideal touchdown point. 

TAFCOS performance is illustrated by the various error measurements of 
figure 5(d). The first six curves show time histories of position, velocity, 
and acceleration errors in the vertical channel. The rough errors are 
between actual and rough- commanded quantities; the smooth errors are between 
actual and smooth commanded quantities. The smooth vertical position error 
indicates how well the aircraft follows the smoothed command, which is the 
actual input to the vehicle control system. It will be seen that the air- 
craft tracking of the smoothed input is extremely good. Maximum smooth posi- 
tion, velocity, and acceleration errors are only 0.15 m (0.5 ft), 0.15 m/sec 
(0.5 ft/sec), and 0.03 g respectively. On the other hand, rough state errors 
are much greater as the aircraft cannot perform the step changes at tip-over. 
These curves emphasize the two problems of trajectory control: first, to 

provide a smooth command signal that is consistent with aircraft capability, 
and second to design a control system that can closely follow the smooth com- 
mand to yield a small error. Figure 5(d) also shows errors in the longitudi- 
nal channel which relate generally to the engine thrust, curve of figure 5(c). 

The curves of figure 5(e) display internal signals to indicate the pro- 
portion of feed-forward and feedback signals that make up the input to the 
trim map. As shown in figure 3, the open-loop feed-forward signal f is 
added to the closed-loop feedback signal f cc to form the specific force- 
command input to the trim map f tc . These signals, which are in velocity axes 
appropriate for entry to the trim map, are shown in the curves of figure 5(e). 
The feed-forward signal f QC has two components: a portion that represents 

the feed-forward acceleration required to follow the flightpath dynamics and 
a portion needed to sustain flight in the presence of gravity. The gravity 
term is the major portion of the feed-forward signal for the vertical axes; it 
should be noted when considering the relative importance of the feed-forward 
and feedback signals. Thus, considering only the variation of the signals it 
is seen that the feed-forward signal is over 90% of the total. The curves of 
figure 5(e) show that for small maneuvers with no atmospheric turbulence, most 
of the control is provided by the feed-forward signals. The differences 
between the trim map and the actual aircraft model show up in the transient 
activity of the first 6 sec of flight where integrators build up to take care 
of the bias effects. 

Figure 5(f) displays second-order effects that were not included in the 
landing run of figures 5(a-e). The A-7E model has almost 1° of hysteresis in 
the throttle actuator, and radar signal information is lost at 150 m (492 ft) 
before touchdown. These effects were observable only in the thrust and in 
some of the longitudinal variables. The touchdown error increased only by 
20 cm (0.66 ft) when these second-order effects were included. 

One other factor that influences touchdown point error is the phase of 
the deck motion at touchdown. Errors are minimum when the touchdown height is 
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near maximum or minimum. Figure 5(g) is a summary plot of the data for 
16 runs of touchdown error as a function of touchdown phase for the condition 
of no disturbances. The deck motion phase is seen to produce a substantial 
touchdown error, because on the 3.5° glide slope an error of 0.3 m (1 ft) 
vertically results in a horizontal error of 5 m (16 ft) . 

Because a distinguishing characteristic of TAFCOS is the calculation of 
the lift and drag coefficients required to follow the trajectory, it is impor- 
tant to compare the commanded coefficients with the values actually achieved. 

A direct comparison of the time histories of the coefficients may be seen in 
figure 5(h). The actual lift coefficient follows the commanded coefficient 
very closely except during the initialization of the first 6 sec and at the 
start of tip-over, when the characteristic aircraft response starts in the 
opposite direction. The drag coefficient likewise follows closely except for 
the abrupt command at tip-over. The variations of commanded lift and drag 
coefficients in figure 5(h) cover only a small part of the trim map that was 
shown in figure 4 and from which it can be seen that the corresponding com- 
manded angle of attack varies from 9° to 14° while commanded thrust coef- 
ficient varies from 0.04 to 0.23. A much greater area of the trim map is 
covered during the gust response runs, as displayed in figures 8(d-f). 


Carrier Approach with Atmospheric Turbulence 

The previous touchdown runs without disturbances are valuable for under- 
standing and analysis of system performance, but in actual service the aircraft 
must operate in an environment of ship-induced and random atmospheric turbu- 
lence. Therefore, a series of 160 simulated landings with disturbances was 
made to collect statistical data on touchdown dispersion. Performance curves 
for a representative run are presented in figures 6(a-f), and a summary of 
touchdown dispersion for the entire series of runs is shown in figure 7 and 
in table 1. To provide a comparison with the previous Navy simulation tests 
at Ames Research Center, the same turbulence and other disturbances were used 
for these TAFCOS simulations. In addition to the usual random wind gusts 
encountered when flying within a few hundred feet of the ocean, a strong pat- 
tern of atmospheric turbulence extends behind the carrier for about 1/2 mile. 
This atmospheric turbulence, shown in figure 6(e), is generated by the bulk of 
the carrier, as it moves through the air at speeds approaching 30 knots. A 
major component of this carrier-wake turbulence or burble is due to the verti- 
cal motion of the carrier in response to wave action. This deck motion, which 
was transmitted to the aircraft for tracking during the final approach, was 
the same as previously shown in figure 5(b). Throttle hysteresis was included 
in all of these runs and its effects can be observed in figure 6(b). Radar 
was inoperative within 150 m (492 ft) of touchdown. 

Most of the quantities shown in figures 6(a-f) were also shown for the 
tip-over run without turbulence in figures 5(a-h). It will be observed that 
the rough and smooth feed-forward commands (6(a)) are nearly the same as 
before. As open- loop commands, they are not dependent on the aircraft 
response to turbulence except through the effect of airspeed. It should be 
recalled that the basic trajectory commands specify position along the 
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trajectory at a commanded airspeed. This effect can be observed by noting 
that the rough- commanded forward velocity relative to the deck in figure 6(a) 
is no longer constant. The effects of disturbances are evident in the veloc- 
ity and acceleration curves of figure 6(a) and in the aircraft variables of 
figure 6(b). Of course, the response to disturbances increases the trajec- 
tory errors. It can be seen from figure 6(c) that the smooth error ordinate 
scales for the vertical channel are two or more times greater than for fig- 
ure 5(d) and that while the errors in figure 5(d) were only at tip-over and 
during deck motion following, they occurred throughout the entire run of 
figure 6(c). 

In the presence of turbulence, lateral response of the aircraft becomes 
significant. As illustrated in figure 6(d), the greatest lateral position and 
velocity responses occur a few seconds before touchdown when strong burble 
effects are encountered. The total specific force command f tc used to com- 
pute the commanded aerodynamic coefficients contains a component from the 
closed-loop command f cc which depends on measured position R and measured 
velocity V; hence, f tc includes the effects of turbulence. The effects of 
turbulence are evident in the commanded thrust coefficient shown in figure 6(b) 
and in the commanded lift and drag coefficients shown in figure 6(f). Consid- 
erably more variation is observed in the actual lift and drag coefficients 
that directly reflect the total turbulence. 

A summary plot of touchdown errors for the entire series of 160 runs is 
shown in figure 7. Runs were made at 16 different values of carrier deck 
motion phase at touchdown — about every 22.5° — as shown by the abscissa of 
figure 7. For each carrier deck motion phase, runs were made with a set of 
10 different turbulence inputs. The same turbulence process operated but with 
a different random number sequence so that roughly the same turbulence varia- 
tions as shown in the top three curves of figure 6(e) were encountered, but 
different peak values of turbulence near touchdown caused a variation of 
touchdown errors. Each point in figure 7 represents the touchdown error for 
one run. The mean error was 3.9 m (+13 ft) (landed long) and the standard 
deviation was 6.7 m (22 ft). This is about half the 12.4-m (40.6-ft) stan- 
dard deviation that resulted from a series of runs in an earlier simulation 
at Ames Research Center in which an operational Navy system was used with the 
same disturbances (ref. 3). It should be noted that on only six runs were the 
landings beyond the 18.3-m (60-ft) limit and none landed more than 18.3 m 
(60 ft) short. There were no ramp strikes as the minimum ramp clearance was 
3 m (10 ft), the mean was 4.9 m (16 ft); and the standard deviation was 0.9 m 
(3 ft). 


Gust Responses 

The two essential functions of any aircraft control system are to maintain 
the commanded trajectory in spite of external disturbances and respond to com- 
mands that modify the trajectory. These capabilities have already been pre- 
sented for the turbulence disturbances and path commands of a standard carrier 
approach. The responses of TAFCOS to much more severe disturbances and to a 
wider range of vertical and lateral commands are presented in the simulation 
results of this section and the next. For the gust responses of this section 
the aircraft has the same system configuration as for the carrier landing. No 
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gains or limits were altered but throttle hysteresis was eliminated. The air- 
craft is commanded to hold an airspeed of 129 knots for a straight and level 
path at an altitude of 152 m (500 ft) . Measurements are made with respect to 
a coordinate system fixed to the Earth. A steady headwind of 30 knots is 
assumed. The gust run for which trajectory variables are displayed in fig- 
ure 8(a) was a series of sharp step gusts. This does not represent a realis- 
tic situation but was chosen to give a concise picture from which the gen- 
eral characteristics of the aircraft response could be obtained. Time-history 
plots of system variables for step gusts of 25 knots sustained for 5 sec from 
six different directions are shown in figures 8(a-c). This simulation was 
accomplished in a single run of 150 sec. 

A 25-knot step gust is a severe disturbance for an aircraft with airspeed 
of only 129 knots; thus, large attitude, flightpath, and thrust corrections 
were required to maintain airspeed and the commanded altitude. However, the 
corrections were generally smooth and nonoscillatory and aircraft angular rates 
and control surface rates were moderate. The throttle and engine power, which 
show characteristic time delays, are very responsive. 

The trajectory variables are shown in figure 8(a) and the aircraft con- 
trols and angular responses are shown in figure 8(b). Time histories of the 
lift and drag coefficients are presented in figure 8(c). Figures 8(d-f) show 
the effect of limiting on the coefficients. 

The first 55 sec of figures 8(a-c) show the response to the down and up 
gusts. There is a somewhat greater trajectory disturbance from the up gust. 
Also the disturbance in forward velocity and airspeed is considerably greater. 
The throttle is driven to its lower limit of 22° and the angle of attack excur- 
sion is greater. 

The tailwind and headwind gust responses, displayed between 55 and 
105 sec, show about the same trajectory disturbances in each direction 
(figs. 8(a-c)). The step in airspeed can be observed as the step gusts are 
applied and removed. It will be observed that the actual lift and drag coef- 
ficients follow the limited commands. This contrasts with the situation for 
the up and down gusts where the actual lift coefficient shows a sharp spike, 
not in response to a command but as a result of the sharp gust which gives an 
instantaneous change of angle of attack. It should be noted that the attitude 
and engine-thrust responses are consistent; that is, they do not act against 
each other by starting in one direction and then reversing one or more times 
during the transient, as can frequently be observed in conventional attitude 
and engine thrust control systems. The tailwind at 55 sec reduces dynamic 
pressure (fig. 8(c)) so lift is lost and altitude drops (fig. 8(a)). An 
increase is then commanded in angle of attack and throttle (fig. 8(b)), so 
that groundspeed (fig. 8(a)) increases and airspeed (fig. 8(a)) recovers 
almost half its commanded value in 5 sec. Vertical velocity (fig. 8(a)) is 
brought almost to zero, and flightpath angle (fig. 8(b)) recovers about half 
its value before the gust abruptly ends at 60 sec. This condition suddenly 
leaves the aircraft at about 15 knots above commanded airspeed, so the 
throttle is quickly reduced to its minimum at 22° and the attitude and air- 
speed reach their normal values before the headwind gust strikes at 80 sec. 
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The response to lateral gusts between 105 and 150 sec in figures 8(a-c) 
shows almost complete symmetry in the two directions. It will be noted that 
there is some horizontal stabilizer response although there had been no aileron 
or rudder response to the longitudinal disturbances. 

The curves of figures 8(d-f) show the effect of trim map perimeter limit- 
ing on the lift and drag coefficients. They are plots of lift against drag 
coefficients for a single 150-sec run to the same scales as the force trim map 
of figure 4 and with the trim map perimeter superimposed. Figure 8(d) shows 
the coefficients C^ c and Cp c that exist before the limiter in figure 10 of 
appendix A. It is seen that because of the severe disturbances of this run, 
values of angle of attack and thrust coefficient are occasionally called for 
that are outside the boundary of the trim map. Figure' 8(e) is a plot of the 
commanded values C-r and C^ c that emerged from the limiter and were sent to 
the trim map. The limits were set to the perimeter of the trim map except 
that in order to provide an angle-of-attack margin, a limit of 16° for a c 
was imposed. It is seen that the curve stays within the limits. Figure 8(f) 
shows the actual coefficients C-^ and Cp that were achieved by the aircraft. 

They follow the limited commands of figure 8(e), except for sharp transients 
when the aircraft responded to gusts. The development of these curves as a 
function of time can be visualized by comparison with the corresponding time 
history plots of figure 8(c). The limiting to desired performance margins is 
quite straightforward for TAFCOS but would be very difficult for a conven- 
tional system. 


Maneuver Command Responses 

The simulation results presented in this section demonstrate the capabil- 
ity of TAFCOS to respond to a complicated set of maneuver commands. They 
roughly correspond to a trajectory that might be required after a carrier land- 
ing wave-off or during a holding-pattern situation. A flightpath of about 
6-min duration over a modified racetrack course, including an altitude change 
of 305 m (1000 ft), was used as shown in figure 9(a). The control system gains 
and limits remained the same as for the carrier landing, but throttle hystere- 
sis was suppressed. Measurements were made with respect to a coordinate system 
fixed to Earth. A steady wind of 30 knots in the negative inertial X-axis 
direction was assumed so that for a constant commanded airspeed of 129 knots 
the ground speed varied from 99 knots to 159 knots. The trajectory consists 
of a series of straight, circular, and helical segments. The input commands 
to the trajectory time sequencer are shown in table 2. They do not quite com- 
mand a closed trajectory, as can be seen by the gap between start and finish 
in figure 9(a). The trajectory consists of 14 segments at a constant commanded 
airspeed. Each line of table 2 gives the flightpath commands for one segment. 
They are converted by the trajectory command sequencer into rough trajectory 
commands. The initial position of each segment is noted by the corresponding 
number along the path in figure 9(a), which is a plot of the trajectory actu- 
ally flown. In an actual flight situation the basic trajectory data could be 
supplied by an air-traffic-control system or entered by the pilot as long as 
the next segment was entered more than 2 sec before the end of the current 
segment. This is necessary as the path segment initialization, discussed in 
appendix A, operates by having the rough command jump to the new trajectory 
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2 sec before the start of the new segment, in the same way that the rough com- 
mands were shown to jump for the tip-over maneuver. 

This is a fundamental characteristic of TAFCOS which allows simple rough 
path segments to be specified with abrupt changes of path angle as in the first 
section of sketch (a) . 


SEGMENT 1 SEGMENT 1 SEGMENT 1 



GEOMETRIC SEGMENT 
JUNCTION 


ROUGH-COMMANDED 

JUNCTION 


SMOOTH-COMMANDED 

JUNCTION 


Sketch (a) 


The trajectory command sequencer then modifies the segment junction by 
introducing a more abrupt discontinuity sooner to provide anticipation, as 
shown in the second section of the sketch (a). This more severe discontinuity, 
but with anticipation, is then smoothed by the filter action of the trajectory 
command generator to give the smooth command shown in the third section of the 
sketch. This final smooth command can be followed quite closely by the air- 
craft and the anticipation gives a smoothly f lared-transition between segments. 

The constant wind makes this trajectory a demanding test; however, it was 
successfully completed, as shown by the performance curves of figures 9(b-d). 
The trajectory variables of figure 9(b) show that the trajectory is smoothly 
followed. There are small disturbances in the vertical velocity, and the air- 
speed shows variations up to 6 knots as turns are made with respect to the 
steady wind. It should be noted that the ground speed varies by ±30 knots when 
flying either directly with or against the wind. 

The variables of figures 9(c and d) show the effects of path changes with 
respect to the wind. The control surfaces have well-damped response, as shown 
in figure 9d. The spoilers are connected to the ailerons but do not deflect 
until the aileron exceeds 3°. The heading angle (fig. 9(c)) is restricted to 
±180°. The jump at about 100 sec does not represent a discontinuity to the 
system since angles are manipulated by direction cosine matrices so all inter- 
nal operations are smooth. 

The rough error curves of figure 9(e) show a substantial sustained error 
between the executed trajectory and the rough commanded trajectory during path 
segment 10. Segment 10 is a downwind turn of 180° between about 200 and 
250 sec. The rough errors reach 213 m (700 ft) for the inertial X axis and 
122 m (400 ft) for the inertial Y axis. For the other segments of the tra- 
jectory, the rough errors are only the usual jump discontinuities at the start 
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of a new segment. The smooth errors are less than 2 m (6.5 ft) in segment 10 
and less than 0.6 m (2 ft) for the other segments, so the aircraft is follow- 
ing the smooth-commanded trajectory quite closely. 

An explanation of the large rough errors in the horizontal plane for seg- 
ment 10 will point out how effectively acceleration limits can be applied to 
the TAFCOS configuration. If the rough- commanded horizontal trajectory i$ 
plotted over the executed trajectory of figure 9(a) the scale is such that no 
difference can be detected. When the rough- commanded and the executed trajec- 
tories are plotted for segment 10 to the enlarged scale of the top curve of 
figure 9(f), however, it can be seen that the executed trajectory does not fol- 
low as tight a turn as the rough- commanded so that the large errors of fig- 
ure 9(e) result. This can be explained by the final three curves of fig- 
ure 9(f), which are time histories of the rough-commanded, smooth- commanded , 
and actual accelerations for the inertial Y axis. The turn of segment 10 
is with a smaller radius (914 m (3,000 ft)) than for the previous turn (1524 m 
(5,000 ft)) and although airspeed is essentially constant over the entire 
trajectory, the 30-knot tailwind gives an inertial velocity of 159 knots at 
the start of the segment-10 turn. These conditions require a radial accelera- 
tion of 7.3 m/sec 2 (24 ft/sec 2 ) which is the peak value of the rough- commanded 
acceleration in figure 9(f). However, internal limits (L3 and L4 of fig. 11) 
were set to allow only 5.8 m/sec 2 (19 ft/sec 2 ) of smooth- commanded accelera- 
tion as shown in the smooth- commanded acceleration curve of figure 9(f). The 
actual acceleration then followed this smooth command very closely. At the 
end of segment 10 the wind has become a headwind and inertial velocity is only 
99 knots. Rough-commanded radial acceleration is only 2.9 m/sec 2 (9.5 ft/sec 2 ) 
so no limiting takes place. Similarly, the ground speed has been reduced to 
129 knots halfway through the turn so the rough-commanded acceleration in the 
inertial X axis is only 4.8 m/sec 2 (15.8 ft/sec 2 ) and no limiting results. 

The limits considered here are limits on commanded acceleration; they are 
discussed in appendix A. (They are not the angle of attack limits mentioned 
in connection with the discussion of the force trim map.) 


PRACTICAL IMPLEMENTATION PROBLEMS 


The simulation of the automatic carrier landing problem using the TAFCOS 
open-loop dynamic trim control concept yielded very promising performance 
results. There are several practical problem areas that must be considered, 
however, before any particular application could realistically be undertaken. 

One of the most important questions concerns the amount of airborne digi- 
tal computer capacity required. A very good estimate of required computer 
capacity for an aircraft installation has recently become available (ref. 2). 
The TAFCOS concept was applied to a commercial aircraft, the deHavilland Twin 
Otter. A successful flight test was completed using essentially the modified 
racetrack course discussed in the section of results. The airborne installa- 
tion used a Sperry 1819A airborne digital computer which has been in flight 
service for 10 years. The TAFCOS installation required only 3000 words of 
memory or about 10% of the 32,000 words available. The airborne computer oper- 
ated at a rate of 20 Hz. The higher frequency rotational dynamics were 
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calculated every cycle and the lower frequency trajectory variables were cal- 
culated every fifth cycle. The TAFCOS calculations used only 10 msec of com- 
puting time per cycle out of an available 50 msec. It was found that the 
required computer capacity was quite reasonable and substantially less than 
that required for a conventional control system with the usual modes of atti- 
tude hold, turn coordination, glide-slope capture, glide-slope tracking, flare, 
and automatic landing. 

The simulation results presented in this report were obtained from an 
IBM 360-67 digital computer and did not simulate the word length or internal 
structure that would be found in an airborne computer; consequently, some per- 
formance differences might be expected due to quantization effects and internal 
computer dynamics. 

Another important question is how to provide the best data for the trim 
map. This simulation relied on wind-tunnel data for the aircraft character- 
istics, but supplementing and correcting those data from flight test results 
would seem desirable. On the other hand, the differences between the trim map 
for this simulation and the aircraft model used for generating the actual state 
variables were largely compensated by integral control. 

A question for the carrier landing is whether to track the deck during the 
final approach or to predict the deck position several seconds ahead and aim 
for that point. A brief study of this problem that showed promise for a pre- 
diction technique was reported in reference 7. The TAFCOS framework is well 
suited to use of a prediction scheme. 

Finally, the question of system reliability is of utmost importance. The 
TAFCOS framework appears to be convenient for reliability analysis and self- 
checking procedures that would utilize its aircraft-response data. 


CONCLUDING REMARKS 


The theoretical basis for a new design concept for computational flight 
control of aircraft has recently been developed at Ames Research Center. This 
report presented the application of the concept to the problem of automatic 
carrier landing. The implementation of the concept was explained by tracing 
the operations in the signal flow path for the system. Simulation results were 
shown for the carrier approach and landing. Further results for extreme dis- 
turbances and a wide range of maneuver cO^nmands were also presented. 

The concept employs a balance of feed-forward and feedback signals. It 
is essentially an aerodynamically based control system that takes full advan- 
tage of a priori knowledge of the aircraft force and moment characteristics 
stored as dynamic trim maps, and utilized by an airborne digital computer to 
calculate the angle of attack and thrust coefficient required to follow a 
trajectory. 
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It was found that the feed-forward signal was about 80% of the total sig- 
nal for most maneuvers and only 20% was needed from the feedback loops. The 
simulation results showed very smooth response to the carrier approach and 
landing commands . The deck motion was closely followed with moderate control 
activity, even for large vertical excursions of the touchdown point. The 
simulation showed that the digital computer requirements could be easily met 
by currently available computers. The TAFCOS configuration allowed straight- 
forward use of limits to control- commanded accelerations and angle of attack 
so that only executable trajectories would be commanded regardless of initial 
path errors, extreme discontinuities at path segment junctions, or severe 
external disturbances. 

It is concluded that the TAFCOS concept potentially has the ability to 
achieve a substantial improvement in landing-point dispersion for the automatic 
carrier-landing task. It also appears that it would be readily adaptable for 
control of the aircraft over a much wider range of trajectory commands. 

Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, October 2, 1979 
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APPENDIX A 
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DETAILED SIGNAL FLOW OF TAFCOS SIMULATION 
NOMENCLATURE 

system matrix — perturbation model 

direction cosine matrix representing rotation from the space axis 
system to the aircraft axis system 

direction cosine matrix representing smooth commanded rotation from 
space axis system to aircraft axis system 

direction cosine matrix representing rough commanded rotation from 
space axis system to aircraft axis system 

estimated value of A 

as 

direction cosine matrix representing rotation from the space axis 
system to the smooth- commanded velocity axis system 

direction cosine matrix representing rotation from space axes to 
rough commanded velocity axis system 

aileron detlection 

commanded aileron deflection 

wing span 

mean aerodynamic chord 
commanded drag coefficient 
commanded lift coefficient 
commanded rolling-moment coefficient 
commanded pitching-moment coefficient 
commanded yawing-moment coefficient 
commanded thrust coefficient 

commanded drag coefficient before limiting 
commanded lift coefficient before limiting 
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h+)> 


C(l) 

C (2) 
C(3) 

M* c ) 
^ 2 ( a c > 
E 2 (r c ) 

E 3(* c ) 

E 3 (-3 c ) 


f c ^^max 
^max 

g 

Gl 9 G2 y • • • 


force coefficient along No. 1 velocity axis 
force coefficient along No. 2 velocity axis 
force coefficient along No. 3 velocity axis 

elementary direction cosine matrix representing a rotation through 
the angle <J> Q about the No, 1 axis 

elementary direction cosine matrix representing a rotation through 
the angle a c about the No. 2 axis 

elementary direction cosine matrix representing a rotation through 
the angle T Q about the No. 2 axis 

elementary direction cosine matrix representing a rotation through 
the angle \p c about the No. 3 axis 

elementary direction cosine matrix representing a rotation through 
the angle -g c about the No. 3 axis 

output of force generation servo 

closed-loop commanded specific force vector 

specific force error vector 

open-loop commanded specific force vector and input to force genera- 
tion servo 

total commanded specific force vector before addition of integral 
term 

total commanded specific force vector 
estimated specific force vector 
time derivative of f c 
second time derivative of f Q 
maximum permissible value of f^(2) 
maximum permissible value of (2) 

• 9 

maximum permissible value of f c (2) 
acceleration of gravity 
gains (three components) 
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I 


identity matrix — three by three 


I 12 identity matrix — twelve by twelve 

J aircraft inertia matrix 

Ll,L2,... limits (three component) 

M aircraft mass 

PLA power lever angle 

PLA C commanded power lever angle 

p estimated aircraft roll rate 

A 23 “ A 32 

Q(A) = A 3 i - A ^3 the column vector function of the rotational matrix A 
A 1 2 “ A 2 1 

q estimated aircraft pitching rate 

q dynamic pressure 

R aircraft position vector 

R c smooth- commanded aircraft position vector 

R e aircraft position vector error between rough and smooth commands 

R e p aircraft position vector error between estimated and smooth command 

R r rough-commanded aircraft position vector 

Rud rudder deflection 

Rud c commanded rudder deflection 

R estimated aircraft position vector 

r estimated aircraft yawing rate 

S wing area 

Splr spoiler deflection 

Splr c commanded spoiler deflection 

Stab stabilizer deflection 
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Stab 

c 

commanded 

stabilizer deflection 



0 

<o(3) 

-03(2) 

S(u>) 

= 

-to (3) 

0 

m(l) the skew symmetric matrix function of the 



to(2) 

-(o(l) 

q angular velocity vector u> 


s Laplace variable 

V aircraft velocity vector 

V c smooth-commanded aircraft velocity vector 

V g aircraft velocity-vector error between rough and smooth commands 

V„ aircraft velocity-vector error between estimated and smooth commands 

V r rough-commanded aircraft velocity vector 

commanded velocity vector with respect to the airmass 

A 

V estimated aircraft velocity vector 

V aircraft acceleration vector 

V smooth- commanded aircraft acceleration vector in space or velocity 

axes as appropriate 

V perturbation acceleration command within trajectory command 

^ generator 

rough-commanded aircraft acceleration vector 

total acceleration command within trajectory command generator 

A 

V estimated aircraft acceleration vector 

A 

W estimated steady wind vector 

a c commanded angle of attack 

commanded sideslip angle 

T c commanded vertical flightpath angle 

F rc rough-commanded vertical flightpath angle 

Af c output of force generation servo — perturbation model 

Clf damping ratio of perturbation model, channel one force servo 

damping ratio of perturbation model, channel one translational 
response 


22 




* 

<i> 

<f> 

i 

$ 

* 


c 

max 

max 

max 

c 

rc 


do 

0) 

(I) 

db 

00 

do 

00 

00 

00 


c 

c 

cc 

f 

oc 

rc 

If 

It 


00 


2f 


00 


2t 


damping ratio of perturbation model, channel two force servo 

damping ratio of perturbation model, channel two translational 
response 

aircraft roll angle 

commanded roll angle 

roll angle limit 

roll angle rate limit 

roll angle acceleration limit 

commanded horizontal flightpath angle 

rough-commanded horizontal flightpath angle 

aircraft angular velocity vector 

estimated aircraft angular velocity vector 

smooth- commanded aircraft angular velocity vector 

total aircraft commanded angular acceleration vector 

closed-loop commanded aircraft angular acceleration vector 

angular frequency in open loop frequency response 

open-loop commanded aircraft angular acceleration vector 

rough-commanded aircraft angular velocity vector 


natural 

frequency of 

perturbation 

model. 

channel 

one 

force servo 

natural 

frequency of 

perturbation 

model , 

channel 

one 

translational 

response 






natural 

frequency of 

perturbation 

model, 

channel 

two 

force servo 

natural 

frequency of 

perturbation 

model. 

channel 

two 

translational 


response 


Subscripts 
c 
e 

ep 


commanded 

error 


P 


perturbation error 
perturbation 
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t 


total 


cc closed-loop commanded 

oc open-loop commanded 

r rough 

rc rough commanded 

tc total commanded 

Other Notations : 

• time derivative 

estimated 

- aerodynamic quantities q, c 

t transpose of an orthogonal matrix, hence also inverse 

solid line denotes 3x1 column vector diagram connection 

= double line denotes 3x3 matrix diagram connection 

INTRODUCTION 


The TAFCOS implementation depends essentially on the airborne digital 
computer. Thus, a complete understanding of all the details that were neces- 
sary for the simulations would require a careful examination of the digital 
computer source programs. Most of the important considerations, however, 
including some that were not covered in the body of this report, can be 
explained by a study of figure 10, which is a diagram of the information flow 
between the major operational elements of the simulation computer program. It 
will be seen that figure 10 has the same basic configuration as figure 3 with 
the addition of several new elements to bring out details. The following dis- 
cussion will trace the information flow in figure 10, commenting primarily on 
the new elements and explaining how the functions of the major elements previ- 
ously described in the body of the report were implemented. 


Aircraft, Control Servos, Sensors and State Estimation 

The aircraft block, the control servo block, the dynamic lag block, and 
the sensor block are not part of the airborne digital computer. The equations 
implemented in the aircraft block represent the complete aircraft as used for 
the simulation of TAFCOS and include a number of secondary aerodynamic effects 
and engine response characteristics that were not taken into account when con- 
structing the trim maps for TAFCOS. The actual aircraft state and trajectory 


24 



variables are measured by various sensors and processed by a state estimation 
system. Data from airborne and shipboard sensors are utilized by such a sys- 
tem. The inputs to TAFCOS are then the hatted quantities from the estimation 
system; they include an estimate of the steady wind velocity vector W and 
the aircraft velocity vector with respect to the carrier V, all expressed in 
the carrier-centered inertial coordinate system defined in appendix B. 


Trajectory Command Sequencer and Segment Initialization 

The trajectory command sequencer contains the desired trajectory geometry 
that describes the aircraft flightpath as a series of straight line, circular, 
and helical segments. This can be prestored information, as in the case of 
the carrier landing where the tip-over the glide- slope geometry is standard- 
ized, or it could be entered a segment at a time, either by the pilot or by 
radio link from an air-traf f ic-control system. The^traj ectory command sequen- 
cer uses the desired airspeed and estimated winds W to calculate an inertial 
acceleration V r which is required to hold airspeed and follow the flightpath. 
This acceleration is integrated to provide commanded velocity V r and posi- 
tion R r . The acceleration command vector is limited to acceptable values, 
thus, although the desired geometric trajectory relative to the carrier is 
commanded, the ground speed and touchdown time may vary due to winds while 
airspeed remains constant. For any desired trajectory, anticipation of the 
next segment is introduced near the end of each segment. Estimated position 
R and velocity V are compared with commanded quantities at the initializa- 
tion of each segment to correct accumulated errors. Therefore, there may be 
small but abrupt changes in command introduced at the start of each trajectory 
segment. The outputs of the trajectory command sequencer are the rough com- 
manded position R r , velocity V r , and acceleration V r , vectors which include 
the carrier deck motion components during the final 12 sec of the trajectory. 


Trajectory Command Generator 

From the preceding discussion of the trajectory command sequencer, it is 
apparent that although the path segments are intended to be flyable the transi- 
tions, possibly discontinuous, are in general not executable. The purpose of 
the trajectory command generator is to generate a smooth, flyable transition 
or capture maneuver. The trajectory command generator employs a four- 
integrator filter on each of three axes to provide smoothing. Information con- 
cerning the proximity of the flight envelope limits is used to ensure the fly- 
ability of the commanded trajectory. 

A conceptual block diagram of the trajectory command generator is shown 
in figure 11. The nomenclature of figures 10-15 generally follows that previ- 
ously employed. Amplitude limits are denoted by the symbol mm /~ ; integrations 
are indicated by the symbol f . The symbol /£ (• ) 2 indicates the square root 
of the sum of the squares of the components of the vector input; that is, its 
magnitude; |V C | also indicates magnitude. The rate limit allows its output 
to change at no more than a specified maximum rate. Figure 11 is really a 
three-channel vector diagram that represents three scalar diagrams. The 
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variables are three-component vectors, such as R r and V c . Rotational trans- 
formations from the inertial space coordinate system to the rough-commanded 

relative velocity axis system are indicated by A vs . It should also be 

rc 

noted that the No. 1 axis of the relative velocity axis system is along the 
direction of the rough commanded relative velocity V r _ft and that the No. 2 
axis is horizontal. The inverse transformation is A vs rc * Gains are indicated 

as Gl, G2, etc. The inputs to the trajectory command generator are rough 
position R r , rough velocity V r , and rough acceleration V r commands, each 
of which has three components. 

The structure of the trajectory command generator is that of a three- 
channel servomechanism in which the plant is a simplified model of the air- 
craft. The state space has dimension 12 and is defined by the outputs of the 
four integrators in three axes. The commanded position of the aircraft is 
defined in inertial coordinates by R c ; the velocity is defined in inertial 
coordinates by V c . The remaining six dimensions are used to model the force 
generation process. The total specific force f c is expressed in the rough 
relative velocity axis system defined by the axis transformation A^ and - 
the final three dimensions define f c . The first channel is identified with 
engine thrust and is primarily associated with the aircraft longitudinal axis. 
The second channel is identified with roll and lateral translational response. 
The third channel is identified with pitch and translational response normal 
to the flightpath. 

The sources of the dynamics in the force generation process are the atti- 
tude and power control systems. A correspondence can be established between 
limits that are to be imposed on engine power, attitude, angular velocity, and 
angular acceleration, and the equivalent limits on the force dynamics. For 
example, a useful approximate link between the lateral force f c (2) and the 
roll angle <j> is simply f c (2 ) = g<j) , where g is the acceleration of gravity. 
This equation, and its derivatives with maximum limiting values of the vari- 
ables, yields the following relations: 

^c^^max ~ S^max 

M^max “ ^max 

■^c(2) max = g^ m ax 

The limiting angular acceleration $ max is related by the moment generating 
process to the aileron and spoiler limits which are given, thus establishing a 
value for f c ^)max - Similarly acceptable maximum values of roll rate <|> 
and roll angle <f> max establish values for f Q (2) x and f c (2) max , respec- 
tively. Limits on the other two axes can be similarly established by asso- 
ciating the No. 1 axis with engine thrust and the No. 3 axis with vertical 
plane maneuvers. 

The limits L5 and L6 in figure 11 are contingency limits on the error 
values R e and V e in space coordinates to cope with unusual cases in which 
extremely large errors could result from faulty commands calling for a trajec- 
tory hundreds or thousands of feet from the present position. They permit 
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such errors to be reduced gradually as the aircraft responds by a prolonged 
transition with acceptable values of various rates and angles. The remaining 
amplitude limits — LI, L2, L3, and L4 — act on quantities in the rough relative 
velocity axis system and can be adjusted along with the gains and rate limit to 
ensure that the desired limits on f Q and its derivatives will not be 
exceeded. 

The force generation process is represented as a simple servo with f Qc 
as input and f c as output and is designed in such a way that regardless of 
the behavior of the driving signal V t the resulting evolution of the force 
is consistent with aircraft limitations. In figure 11, rough- and smooth- 
commanded position R c , rough- and smooth-commanded velocity V r , V , and 

rough- and smooth- commanded acceleration V r , V c are in inertial coordinates; 
the force variables are in the rough relative velocity axis system. The rough 
acceleration command V r has the perturbation component Vp added to give 
the total acceleration command V t . The addition of a force command to over- 
come gravity yields the specific force f as input to the force generating 
servo. In the steady state, f c will equal f QC because of the unity negative 
feedback. The gains G3 and G4 are chosen to give a desired natural frequency 
and damping ratio to the dynamic buildup of f c . This commanded specific 
force f Q is then resolved to inertial axes and the counter gravity force 
removed to give the total commanded path acceleration V Q in inertial coordi- 
nates. Simple integration then gives velocity V c and one more integration 
gives position R c . 

It might be expected that the specific force f c would serve as the input 
to the trim map; however, f QC is used as the open-loop specific force command 
because it provides the required dynamic anticipation of f Q while being equal 
to it in the steady state. 

Open-loop excitation of the trajectory command generator is # provided by 
the acceleration command from the trajectory command sequencer V r . In addi- 
tion, regulation feedback Vp is provided to force the trajectory command 
generator to follow the trajectory from the sequencer asymptotically. Several 
limits are included, as shown in figure 11, to ensure stable, well-behaved 
regulation for large initial errors. 

The trajectory command generator provides the open- loop specific force 
command f QC , which is combined with the closed-loop specific force command 
from the trajectory regulator and sent to the force trim map section further 
down the hierarchy of the control logic, where it is converted to correspond- 
ing attitude and power commands. Additional outputs of the trajectory com- 
mand generator are R c and V c ; they describe the reference trajectory and are 
sent to the trajectory perturbation regulator for comparison with the corre- 
sponding measured quantities and subsequent regulator action to reduce trajec- 
tory errors. Since TAFCOS is primarily a feed-forward system that relies only 
partially on feedback, the regulator is properly a perturbation regulator. 

The following paragraphs provide a more detailed explanation of how the 
gains in figure 11 were chosen by using linear control theory. Because fig- 
ure 11 is a vector block diagram that represents three channels, the variables 
are three-component vectors; that is: 
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R r 


R 


XI 
R r2 
R r3 

and the gains were chosen to be three by three diagonal matrices; that is: 

i 0 0 

G\ ~ 0 ^12 0 

0 0 Gj 3 


The direction cosine matrix in figure 11 represents a rotation from 

inertial axes to the rough- commanded relative velocity axes system whose first 
axis is lined up with the direction of the rough- commanded relative-velocity 
V r - W and whose second axis is horizontal. This matrix differs slightly 
from A vs in figure 13 which represents a rotation from inertial axes to the 
smooth- commanded relative-velocity axis system whose first axis is lined up 
with the direction of the smooth- commanded relative velocity V c - W and whose 
second axis is horizontal. This distinction is necessary to avoid an internal 
feedback loop (fig. 11) that would arise if V Q , which is computed in the loop, 
were used instead of V r . 

The perturbation model of the command generator for a linear analysis is 
shown in sketch (b). The symbol Af c is introduced instead of V c because 
the addition and subtraction of the constant gravity term has been ignored. 


FORCE GENERATION SERVO 



Sketch (b) 
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The integrator outputs are chosen as the state variables of this twelve- 
dimensional system. The unforced state equations — in vector matrix form — 
are: 


R 


0 10 0 


R 

c 




c 

V 


0 0 10 


V- 

c 

_ 



c 

A ^c 


0 0 0 1 


Af c 



” G 1 G 3 ~ G 2 G 3 ~ G 3 _G 4 G 3 


f 

c 


where I is the three-by-three identity matrix and the other quantities are 
vectors and matrices as previously defined. The characteristic equation of the 
system is det(sl 12 - A) = 0 , where A is the system matrix and I 12 is the 
twelve-by-twelve identity matrix. This leads to the determinant equation 


which yields 


Is 

-I 

0 

0 

0 

Is 

-I 

0 

0 

0 

Is 

-I 

G3G1 

G 3 G 2 

g 3 

(Is + G3G4) 


Is 4 + G3G4S 3 + G3S 2 + G 3 G 2 S + G3G! = 0 


This vector equation represents a three-axis decoupled system because the gain 
matrices are diagonal. Thus, each of the corresponding scalar equations rep- 
resents a single channel. 

The roots of the scalar characteristic equation are the poles of the 
transfer function for the corresponding channel. The gains for the various 
channels in sketch (b) can be determined by any of the conventional pole place- 
ment techniques. The following discussion explains in detail how they were 
determined for the second channel. The force generation servo represented 
essentially by the first two integrators was specified to have a second-order 
response with natural frequency = 2.85 rad/sec and damping ratio 

^ 2 f = The lateral translational response generally associated with the 

last two integrators was specified to have a second-order response with natu- 
ral frequency 0J 2t = 0.67 rad/sec and damping ratio £ 2t = 0 . 71 . This leads 
to a desired transfer function that is a cascade of two second-order systems. 
The denominator of this transfer function then determines the characteristic 
equation: 


(s 2 + 2£ 2f 0) 2f s + 0) 2 f ) (s 2 + 2? 2t 0J 2t s + w 2 t ) 

When this desired characteristic equation is set equal to 
equation of sketch (b) we have: 


- 0 

the characteristic 
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S 4 + s3(2 ^2f“2f + 2e 2t U, 2t ) + s2(45 2f 5 2t“2f“2t + *2f + U >2t ) 


+ s(2 ?2f (o 2fl u| t + 25 2 t « 2 t a.| f ) + = s* + s 3 G G + s^G, + sG 3 G 2 + G 3Gl 


A direct matching of the coefficients gives the gains G 1 = 0.29 , G^ = 0.76, 

G 3 = 12 . 6 , G 4 = 0.41. 

The same gains were also used for the third channel. For the first chan- 
nel, which is identified with the slower engine thrust response, the desired 
parameters were taken as a) lf = 1.52, = 0.75, u>, t = 0.53, ? lt = 0.73; these 

parameters yield the gains: Gj = 0.15, G 2 = 0.56, 63 =4.9, and G 4 = 0.7. It 

is apparent that all of the techniques of modern control theory and optimiza- 
tion could be used at this point in the trajectory command generator design 
to give specified dynamics for each channel. One interesting feature of 
sketch (b) is the inclusion of the gain G 3 in the direct path. This pre- 
vents the representation from being exactly the controller canonical form, but 
permits unity feedback to close the force generator servo loop. This gives a 
direct physical representation for convenient numerical checking during system 
tests and observation of limiting behavior. 

The design described above was examined by classical servomechanism analy- 
sis. The loop was opened at point P in sketch (b) to give the servo config- 
uration of sketch (c) which has an open-loop transfer function of 

Q = G 3 (Gj 4- G 2 S) 

P s 2 (s 2 4- G 3 G 4 s + G 3 ) 



Sketch (c) 
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This open-loop transfer function, of course, has different roots and poles than 
the specified closed loop. The amplitude ratios and phases of the open-loop 
frequency response and of its factors are shown in figure 12 plotted against 
the logarithm of angular frequency oj^. Curve A is the amplitude ratio of the 
force servo factor 


F = ^3 = 1 

P s 2 + G 3 G 4 S + G 3 1 + G^s + (s 2 / G 3 ) 


with break point at a)f = 3.55. In figure 12, curve B is the corresponding 
phase; curve C is the amplitude ratio of the remaining factor 




■) 


with break point at u)£ = 0.377. The corresponding phase is curve D. The 
total open-loop amplitude ratio from P to Q is curve E and the corresponding 
phase is curve F. It can be seen that this yields a very satisfactory Bode 
plot with a phase margin of 46° and a gain margin of 16 dB. 


Time domain responses were then investigated on the simulation and iter- 
ated to give the most satisfactory transient performance. The final gains thus 
selected are those used for the example calculations; this accounts for the 
three-figure precision. 


Trajectory Perturbation Regulator 

The trajectory command generator discussed in the preceding section pro- 
vides an open-loop control and the corresponding reference trajectory. Because 
the aircraft equations of motion and the environment are never represented in 
every detail, a closed-loop feedback is needed to regulate the errors caused 
by deviation of the real process from the model used in the command generator. 
The purpose of the trajectory regulator is to provide control of the trajectory 
errors. The structure of the regulator is shown in figure 13. The structure 
of figure 13 is quite similar to the first section of figure 11. Note that the 
position and velocity error signals are formed by a comparison of smooth com- 
mands R c and V c , from the trajectory command generator, with measured values 
R and V, rather than by comparison of rough and smooth commands. Figure 13 
also includes an integrator section that was not shown in figure 10. The tra- 
jectory errors in position, R ep = -R 4- R c , and velocity, V ep = -V + V c , which 
are in inertial coordinates, are transformed into velocity axes, passed through 
limits and gains, and then combined to form a closed-loop corrective perturba- 
tion specific force command f cc which is added to the open-loop command f QC 
from the trajectory command generator. The gains and limits, applied in the 
velocity axis system, are designed by assuming the same model of the force 
generation process that was used in the trajectory command generator. In addi- 
tion to the position plus (velocity) feedback, integral feedback is used as 
shown in the figure. This loop is applied to the force generation process, 
and it unloads the corrective acceleration f cc . The effect of integral con- 
trol of the error (f £ = f pc - f) is limited by L9 to 0.1 g. The estimated 
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specific force f is obtained from the accelerometer package. Consequently, 
slow errors in the a priori model of the force generation process resulting in 
acceleration errors up to 0.1 g will cause no trajectory offsets. The result- 
ing specific force f tc representing open-loop command plus closed-loop per- 
turbation is next expressed in coefficient form. 


Force Coefficient Development 

The total commanded specific force vector f tc in figure 10 is in rela- 
tive velocity axes. It is multiplied by M/qS to provide force coefficients 
C(l), C(2), and C(3). The drag acts along the negative_No. 1 a^is in this sys- 
tem so the commanded drag coefficient before limiting Cd c is merely the neg- 
ative of the No. 1 force coefficient. The total lift coefficient before limit- 
ing as defined in re ference 2 is d ue to all forces perpendicular to the drag. 
Therefore is /c(2) 2 + C(3) 2 • Also, the commanded roll angle <j> c is 

calculated as the arc tangent of the ratio of the second and third force coef- 
ficients. At this point in figure 10 the commanded lift and drag coefficients 
before limiting C L and Cp have been generated in response to the open-loop 
commands supplemented by the closed-loop commands so that for severe atmo- 
spheric disturbances they may be beyond the aircraft capability. Therefore, at 
this point they are conveniently limited to the perimeter of the force trim 
map of figure 4, or to some sub-envelope based on desired stall margin, to 
give the limited commanded coefficients and which are then sent to 

the trim map. The operation of these limits can be c seen dramatically in 
figure 8(e) compared to figure 8(d). 


Force Trim Map and Thrust Calculations 

The aircraft model for this simulation represents only the landing config- 
uration; hence, no effects of changing flaps are included. Thus, the simple 
two-dimensional trim map of figure 4 could be developed. Furthermore, the trim 
map includes only the major aerodynamic forces due to the angle of attack and 
stabilizer. Other secondary aerodynamic lift, drag, and side forces, used in 
the complete aircraft model, are neglected for the force trim map but could be 
included if more complete open-loop control fidelity was desired. Some of 
them are in fact included in the moment trim map calculations. As noted in 
the body of the report, a three-dimensional trim map to include flap angle con- 
figuration changes was developed for another application. 

The commanded thrust coefficient Ct c from the force trim map is multi- 
plied by qS to give the total required thrust. The digital computer then 
interpolates a table of thrust vs power-lever angle to determine the required 
throttle setting. This commanded power-lever angle PLA C is then modified by 
a dynamic lag and hysteresis to simulate the action of the throttle actuator 
that is used in the A-7E aircraft. 

In a different application, in which engine response time and different 
nonlinear processes that occur during thrust buildup and decay were particu- 
larly important, a thrust perturbation regulator was developed that used 
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measured values of engine speed feedback to provide a closed-loop perturbation 
command similar to the action of the trajectory perturbation and the attitude 
perturbation regulators. 


Angle and Angular Rate Calculations 


As explained in appendix B, the rough- commanded aircraft attitude matrix 
Aas rc is constructed as the product of five elementary direction cosine 

matrices 


‘as 


rc 


E 2 (a c )E 3 (-B c )E 1 ((j) c )E 2 (r Q )E3(if; c ) 


The commanded flightpath angles and T c are computed from the components 

of the commanded^velocity with respect to the steady wind. The estimated 
steady wind is W and the commanded inertial velocity is V c so that the 
vector relative velocity is V rw = V c - W and the commanded flightpath angles 
are computed as: 


ib 

Y c 


tan 1 


-Vrw (2) 

W 1 * 


r c = tan 


-1 


- v rw< 3 > 


/v (l) 2 + V (2) 2 
rw v 7 rw v 7 


If desired for a particular trajectory, a commanded angle of sideslip B c 
can be generated by the trajectory command sequencer. The rough commanded 
angular rate vector co rc is shown in appendix B to be: 


CQ 


rc 


-V c (2) 


tan r c 


-V c (3) 

l V cl 

V c (2) 

l V d 


where V c is the smooth- commanded acceleration vector expressed in velocity 
axes for this particular calculation and | V c | is the magnitude of the smooth- 
commanded velocity vector. 


Attitude Command Generator and Perturbation Regulator 

The variables in the attitude command generator and perturbation regula- 
tor include vectors and matrices. The details of the implementation can be 
seen from figure 14. They are manipulated and integrated in much the same 
way as described in the detailed discussion of figure 15 in appendix B. 
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Commanded Moment Coefficients and Moment Trim Map 

The commanded moment vector is first obtained from the commanded angular 
acceleration m c by matrix multiplication by the inertial matrix J. Then the 
moment coefficients are obtained by dividing the pitching torque component by 
qSc and the roll and yaw torque components by qsb. The moment trim map 
involves both tabular interpolations and algebraic equation solutions. For the 
lateral response, the roll and yaw torque contributions due to commanded angle 
of attack, roll, and yaw rates are calculated from tabular data. The remaining 
portion of the torques required to produce the commanded roll and yaw accelera- 
tions can be expressed in two simultaneous linear algebraic equations in the 
lateral controls; the equations are solved for the required rudder, and lateral 
stick deflection. The aileron and spoiler are then given as functions of the 
lateral stick. The pitching torque contribution due to the commanded angle of 
attack, aileron and spoiler deflections, and commanded pitch rate is then cal- 
culated using tabular data. It is subtracted from the total required pitching 
torque to find the torque contribution needed from the horizontal stabilizer. 

A table of torque vs commanded angle of attack for various stabilizer deflec- 
tions is then interpolated to find the required stabilizer deflection. 

This completes the description of the calculations carried out for the 
simulation of TAFCOS. These calculations must be carried out by the airborne 
digital computer 20 times per second. To provide insight into the system, and 
for convenience and flexibility in the simulation, a number of special subrou- 
tines were employed to execute vector and matrix computations and for one- and 
two-dimensional table interpolations as well as for recording computer output. 
However, when the TAFCOS concept was flight tested, the airborne computer 
requirements were reduced to a minimum. Matrix calculation subroutines were 
replaced by efficient in-line coding, tabular data were replaced by piecewise 
linear analytic expressions, and slow and fast loops were employed. 
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APPENDIX B 


MATHEMATICAL DETAIL 
NOMENCLATURE 


A = 


A 


as 


a l 1 

a 12 

a l 3 



a 2 1 

a 22 

CO 

CM 

cd 

direction 

tional 

cosine matrix representation of a rota- 
transformation 

a 3 1 

a 32 

a 3 3 




direction cosine matrix representing rotation from space to body axes 


A as direction cosine matrix representing the rotation from space to rough 

commanded aircraft body axes 


A 


vs 



as 


dA 

dt 

= A 


direction cosine matrix representing the rotation from space to veloc- 
ity axes 

time derivative of the matrix A 

s inverse and transpose of the orthogonal matrix A^ s 


a ll> cosine of angle between axes of two coordinate systems as indicated by 

a 32 , the subscripts 

etc. 


E x (y) elementary direction cosine matrix representing a rotation through the 
angle y about the axis x. For example: 


E 2 (ct) 


cos a 0 -sin a 

0 10 

sin a 0 cos a 


h a angular momentum vector of the aircraft in body axes represented as: 

h (1) 

V« 

h (3) 
a v ' 


h g angular momentum vector of the aircraft in space axes 

I the identity matrix 3x3 

matrix representation of the inertia operator of the aircraft in body 
axes 3x3 
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J- 1 

a 


M 


M 


S(w) = 


V. 


r 

|V. 


3 

r 

r c 

$ 

V 

* 

ib 

y c 

0) 


inverse of the matrix 

a 

matrix representation of the inertia operator of the aircraft in space 
axes 

total torque vector acting on the aircraft in body axes 
total torque vector acting on the aircraft in space axes 
0 oo(3) —oo (2) 


—oo (3) 0 oo(l) 

oo (2 ) -0)(1) o 


the skew symmetric matrix function of the 
vector oo 


smooth-commanded velocity vector in velocity axes 
rough- commanded velocity vector in space axes 
smooth- commanded acceleration vector in velocity axes 
rough-commanded acceleration vector in space axes 
magnitude of V c 
angle of attack 
commanded angle of attack 
sideslip angle 
commanded sideslip angle 
vertical flightpath angle 
commanded vertical flightpath angle 
Euler angle about No. 1 axis 
roll angle 

commanded roll angle 
Euler angle about No. 3 axis 
horizontal flightpath angle 
commanded horizontal flightpath angle 
angular velocity vector associated with the rotational matrix A rep- 

UiV 


resented as: 


to(2) 

w(3) 


36 



0 


Euler angle about the No. 2 axis 


co a angular velocity vector of aircraft axis system (associated with A as ) 

expressed in aircraft body axis system 

o) c smooth-commanded angular velocity vector (associated with A ag ) in 

body axes c 

°°rc rough-commanded angular velocity vector (associated with A^) in 

velocity axes 

Subscripts 


a 

aircraft body axis system 

c 

commanded 

r 

rough 

s 

space axis system 

V 

velocity axis system 

Other 

Notations 

• 

time derivative 

-1 

inverse of a matrix 


COORDINATE SYSTEMS 


TAFCOS requires the definition of several coordinate systems. All are 
right-handed, three-axis orthogonal systems. Vectors are expressed by their 
three components in the appropriate system. Sketch (d) illustrates the chain 
of coordinate systems that transforms quantities from inertial space coordi- 
nates to aircraft body axes. The inertial space axis system is centered at the 


AXIS 


INERTIAL 

SPACE 


HEADING 


VELOCITY STABILITY 


WIND TUNNEL AIRCRAFT 

BODY 



Sketch (d) 
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carrier’s undisturbed center of pitch with its first axis along the carrier 
deck and its third axis downward. For the short times involved in the landing, 
the carrier velocity has been assumed to be constant over a flat nonrotating 
Earth with the steady wind along the deck; hence, this is an inertial system. 
The aircraft velocity vector with respect to relative wind is along the first 
axis of the velocity axis system. The velocity axis system thus results from 
a rotation of the space axis system about its third axis through the horizontal 
flightpath angle \jj followed by a rotation about its second axis through the 
vertical flightpath angle T. The next rotation about the No. 1 velocity axis 
through the roll angle <j> defines what has been called the wind-tunnel axis 
system. Rotation about its third axis through the sideslip angle (negative 3) 
leads to the stability axis system, and a final rotation about its second axis 
through the angle of attack a defines the aircraft body axis system. It 
should be noted that this series of five rotations achieves the same aircraft 
body axis orientation as the three conventional rotations through the Euler 
angles 0, and Three of these axis systems are used extensively. The 

space axis system (s) is used for basic trajectory definition and for differen- 
tiation when inertial accelerations are involved. The velocity axis system 
(v) is particularly important because the aircraft center of gravity is fixed 
at the origin and the relative velocity vector is aligned along the No. 1 axis. 
The aircraft lift and drag are easily expressed in this system. The aircraft 
body axis system (a) is used for developing forces and torques due to the 
engine and control surfaces and for body rate calculations. The coordinates 
of a vector are transformed from one of these coordinate systems to another by 
multiplication by the direction cosine matrix between the systems. The first 
subscript of the rotation matrix symbol A represents the final system, and 
the second subscript represents the initial system. For example A^g is the 
rotation matrix from space to velocity axes. The rotation from space to air- 
craft body axes is A as = E 2 (a)E 3 (-3)E 1 (4>)E 2 (T)E 3 (i/0 . It should be noted that 
this does not represent the rotation of a vector. The vector remains fixed 
but is represented in various coordinate systems which are rotated with 
respect to each other. 

The Skew Symmetric Matrix Function of a Vector 


The vector cross product a) x v where oj and v are 


tors is equivalent to the matrix product 
metric matrix 

-S (oj) v 


0 

co 3 

— co 2 

S(u)) = 

-(0 3 

0 

COl 


co 2 . 

-0)2 

0 


three-component vec- 
S(o)) is the skew sym- 


This can be shown by forming the vector cross product and comparing term-by- 
term with the matrix product. The skew symmetric matrix function S(w) is in 
some sense the matrix equivalent of the vector m and is used in relations 
where both three component vectors and three-by- three matrices are involved. 
The skew symmetric matrix function of the angular velocity vector appears in 
figures 14 and 15. 
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The Vector Function of a Matrix 
Let the rotation matrix A be 


A = 


a ll 

a l 2 

a l 3 

a 2 1 

a 22 

a 2 3 

a 3 1 

a 32 

a 3 3 

is defined 

as 

x i 

a 2 3 

a 32 

2 

a 3 1 

“ a l 3 


a l 2 

a 21 


In some sense Q(A) is the vector equivalent of the matrix A as can be 
seen by noting that under appropriate conditions the Q and S functions can- 
cel. In general Q ( S (cn) ) = w and for a skew symmetric matrix A, S(Q(A)) = A 
The vector function of a rotational matrix Q(A) appears in figure 14 where it 
is used to convert a matrix representing an angular error into a vector that 
can be added to another vector representing an angular velocity error to form 
a total error vector. 


Angular Rate Calculations 

The time derivative of a direction cosine matrix A 
where S(u)) is the skew symmetric matrix 


is: A = S(m)A 


S(w) = 


0 

-03(3) 

o)(2) 


o) (3) 
0 

-03(1) 


-03(2) 

03(1) 

0 


whose elements are the components of the angular velocity vector associated 
with the rotation, that is: i 

03 ( 1 ) 

03 = ! U) ( 2 ) 

03(3) 

The rough- commanded aircraft orientation matrix (A aSrc ) is the input to 
the attitude command generator in figures 2, 3, and 10. It is desirable to 
supplement this open-loop commanded attitude by the corresponding open-loop 
commanded attitude rate. 

This would strictly be the angular velocity vector associated with A as 
This in turn depends on the time derivatives of all five angles involved in 


rc 
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the transformation. However, only the major components were implemented for 
this application. Hence, only the time derivatives of ip c and F c were used 
to command open-loop attitude rate. The commanded flightpath angles ip c and 
r c , which define the orientation of the relative velocity axis system, may 
change over a wide range during an extended period of time as the flightpath 
is described. Thus, their time derivatives can maintain substantial, slowly 
changing values. Furthermore, these commanded angles are computed from the 
components of the commanded aircraft velocity with respect to steady wind and 
hence are noise free. On the other hand, commanded angle of attack a Q 
changes rapidly to initiate maneuvers, contains a noise component due to feed- 
back from the trajectory regulator, and has zero average rate of change over 
relatively short periods. Commanded sideslip 0 C , which would be needed for 
crabbing and decrabbing, was not used in these simulations and in any event 
would have only short duration, nonzero values of time derivative. Commanded 
roll angle cj> c about the velocity axis likewise is characterized by rapid 
changes of short duration. The angular velocity vector of the velocity axis 
system is given in terms of the angle rates ip c and F c as: 

~h sin r c 

“rc = r c 

\b cos r n 
c c 

Now \p c and F c are the flightpath angles generated to command the air- 
craft center of gravity to follow the trajectory. The trajectory has been 
defined in terms of the corresponding accelerations. The rough commanded 
acceleration V r (fig. 10) is expressed in space axes at the input to the tra- 
jectory command generator and produces the smooth commanded acceleration V c 
in velocity axes. The flightpath angle rates and ip c depend on components 

of this acceleration that are perpendicular to the velocity vector and on the 
magnitude of the velocity |V C | : 

. _ V c (3) _ V c (2) 

c Kl ^ I v c i cos r c 

Substituting these values into the expression for rnrc gives the rough 
commanded open-loop angular velocity, which is used as an input to the attitude 
command generator in figure io. 


The Matrix Form of the Rotational Dynamic Equation 

The rotational dynamic portion of the aircraft simulation represented in 
the aircraft block at the extreme right end of figure 3 uses a matrix formula- 
tion of the dynamic equation derived below. Let M g and M a be the total 
torque — aerodynamic and propulsive — expressed in the inertial frame and the 
body frame respectively. Let h g and h a be the angular momentum vector in 
space and body coordinates. Let J g and J a be the inertia matrix in space 
and body axes. Let cu a be the angular velocity vector of the body axes with 
respect to the space axes expressed in body axes coordinates and let A as 
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be the rotation matrix from the space axes to the body axes. Note that the 
vector quantities are all three-element column vectors; that is. 


^a 



and the matrices are three-by-three arrays. Now 


= 


= 


A as M s 

A as^s 


as 


= S(w a )A 


as 


from a preceding section of this appendix. Also by neglecting spinning rotors 
we have 


h a = J a 03 a 

and because J a is constant in the body frame 

\ = J a*a 

Furthermore, the basic dynamic relation in inertial axes is 

h = 


The relation using the easily available quantities M a and h a may now be found 
by the following manipulations: 


then 


^a ^as^s 


h a = A as h s + A as h s = J a“a 


“a = ^'^^a^as^a + A as M sJ 


i a = J I 1 [ S (“a )h a + M a J = V [M a + s (“ a ) J a w a] 


A block diagram of this expression is included in the first section of 
figure 15. 


The final portion of figure 15 illustrates how the angular velocity vec- 
tor w a is integrated to obtain the aircraft attitude expressed as the direc- 
tion cosine matrix A ag . The vector oo a is first used to construct the skew 
symmetric matrix function S (o> ) . Then the relation for the derivative of a 
direction cosine matrix A ag = s (^ a ) A as is implemented by the matrix 
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multiplication indicated. This is followed by a straightforward integration 
to give A as . The Gram-Schmidt orthogonalization is then required, otherwise 
cumulative computer round-off error would soon lead to a nonorthogonal matrix 
in which individual elements exceed unity. Angles are treated throughout the 
system by manipulation of their direction cosine matrices. This completes the 
aircraft rotational dynamics. However, for purposes of angular readout during 
the simulation, a computer routine has been developed (ref. 8) that further dis- 
plays the Euler angles corresponding to the direction cosine matrix A 00 : 

A as = E 1 ($)E 2 (0)E 3 (T) 


42 



REFERENCES 


1. Meyer, G. ; and Cicolani, L. S.: A Formal Structure for Advanced Automatic 

Flight-Control Systems. NASA TN D-7940, 1975. 

2. Wehrend, W. ; and Meyer, G. : DHC-6 Flight Tests of the Total Automatic 

Flight Control System (TAFCOS) Concept. NASA TP-1513. (In preparation) 

3. Smith, G. Allan; McNeill, Walter E.; and Huff, Richard B.: A Comparison 

of Restructured Designs for the Automatic Carrier Landing System. Pro- 
ceedings of Flight Control Systems Criteria Symposium, U.S. Naval Post 
Graduate School, July 1978. 

4. Smith, G. Allan; and Meyer, George: Total Aircraft Flight-Control 

System — Balanced Open- and Closed-Loop Control With Dynamic Trim Maps. 
Proceedings of the Third Digital Avionics Systems Conference, Ft. Worth, 
Texas, Nov. 1979, IEEE Catalog No. 79CH1518-0. 

5. Cicolani, L. S.; Sridhar, B.; and Meyer, George: Configuration Manage- 

ment and Automatic Control of an Augmentor Wing Aircraft with Vectored 
Thrust. NASA TP-1222, 1979. 

6. Meyer, G.: On the Use of Euler’s Theorem on Rotations for the Synthesis 

of Attitude Control Systems. NASA TN D-3643, 1966. 

7. Sidar, M. ; and Doolin, B. F. : On the Feasibility of Real-Time Predictions 

of Aircraft Carrier Motion at Sea. NASA TM X-62454, 1975. 

8. Meyer, G. ; Lee, H. Q. ; and Wehrend, W. R. , Jr.: A Method for Expanding a 

Direction Cosine Matrix into an Euler Sequence of Rotations. NASA 
TM X-1384 , 1967. 


43 



TABLE 1.- SUMMARY OF 160 SIMULATED APPROACHES WITH TURBULENCE. 


Condition 

Dispersion 

Maximum deck displacement due to pitch and heave 

2.4m (8ft) 

Mean touchdown error (long) 

3.9 m (13 ft) 

Standard deviation of touchdown error 

6.7 m (22 ft) 

Mean ramp clearance 

4.9 m (16 ft) 

Standard deviation of ramp clearance 

0.9 m (3 ft) 

Touchdowns beyond allowable area 

6 

Touchdowns short of allowable area 

0 


TABLE 2.- CLOSED TRAJECTORY PARAMETERS. 


Segment 

number 

Trajectory type 

1 

Straight level 

2 

Straight climb 

3 

Helical 


4 

Straight 


5 

Helical \ 

f 

6 

Circular level 

7 

Straight 


8 

Circular 


9 

Straight 


10 

Circular 


11 

Straight 1 

f 

12 

I descend 

13 

descend 

14 

j level 


Climb 

angle, 

deg 

Path length 

m 

ft 

— 



512 

1680 



349 

1146 







914 

3000 



1755 

5760 



3947 

12950 

\ 

' 

628 

2060 

-6 


2509 

8230 

- 0 . 

9 

1161 

3810 

— 

— 

610 

2000 


Duration. 


Radius 


deg 

m 

ft 

90 

1524 

5000 

64 

1524 

5000 

71 

1524 

5000 

45 

1524 

5000 

180 

914 

3000 
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Figure 1.- Aircraft carrier landing situation. 
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Figure 2.- 


Basic open-loop feed-forward control concept. 



Figure 3.- Total aircraft flight-control system (TAFCOS) 
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Figure 4.- Force trim map. 
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(a) Trajectory variables. 

- Tip-over run without turbulence. 
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(a) Trajectory variables - Concluded. 


Figure 5.- Tip-over run without turbulence - Continued. 
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(b) Vertical deck motion variables. 

Figure 5.- Tip-over run without turbulence — Continued 
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(c) Aircraft variables. 

Figure 5.- Tip-over run without turbulence - Continued. 
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(d) Trajectory errors. 

Figure 5.- Tip-over run without turbulence - Continued. 
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(e) Internal feed-forward and feedback signals. 
Figure 5.- Tip-over run without turbulence - Continued. 
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(f) Radar signal lost at 150 m (493 ft) before touchdown, throttle hysteresis 
included. Plots of major variables that were noticeably affected. 

Figure 5.- Tip-over run without turbulence - Continued. 
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(g) Touchdown error as a function of deck-motion phase. 
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(h) Commanded and actual lift and drag coefficients. 
Figure 5.- Tip-over run without turbulence - Concluded. 
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(a) Trajectory variables. 


Figure 6.- Tip-over run with turbulence. 
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(a) Trajectory variables - Concluded. 


Figure 6.- Tip-over run with turbulence - Continued. 
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(b) Aircraft variables. 

Figure 6.- Tip-over run with turbulence - Continued. 
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(c) Trajectory errors. 

Figure 6.- Tip-over run with turbulence - Continued. 
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(d) Lateral variables. 

Figure 6.- Tip-over run with turbulence - Continued. 
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(e) Turbulence and radar noise inputs. 
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(f) Commanded and actual lift and drag coefficients. 


Figure 6.- Tip-over run with turbulence - Concluded. 
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Figure 7.- Landing position errors for 160 runs with turbulence. 
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(a) Trajectory variables. 

Figure 8.- Level run with step gust disturbances. 
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(a) Trajectory variables - Concluded. 


Figure 8.- Level run with step gust disturbances - Continued 



HORIZONTAL 
STABILIZER, deg 


POWER-LEVER 
ANGLE, deg 


THRUST, lb 


PITCH ANGLE, deg 




(b) Aircraft variables . 

Figure 8.- Level run with step gust disturbances - Continued. 
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(b) Aircraft variables - Concluded. 

Figure 8.- Level run with step gust disturbances - Continued. 
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(c) Lift and drag coefficients. 

Figure 8.- Level run with step gust disturbances - Continued 
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(d) Commanded lift vs drag coefficients, not limited. 
Figure 8.- Level run with step gust disturbances - Continued 
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(f) Actual lift vs drag coefficients. 
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(a) Flight profiles. 

Figure 9.™ Closed-course trajectory. 
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(b) Trajectory variables. 
Closed-course trajectory - Continued. 
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(c) Angle variables. 

Figure 9.- Closed-course trajectory - Continued. 
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(d) Aircraft variables. 

Figure 9.- Closed-course trajectory - Continued. 
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(e) Position errors. 

Figure 9.- Closed-course trajectory - Continued 
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(f) Downwind turn details. 

Figure 9.- Closed-course trajectory - Concluded. 
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Figure 10.- TAFCOS simulation information flow. 
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Figure 12.- Open-loop frequency response of traj ectory-command-generator 

perturbation model. 
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Figure 15.- Rotational dynamics in aircraft section of simulation 
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